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ABSTRACT 


The shoreward and seaward propagating tidal wave eomponents were determined using 
four co-loeated pressure and velocity sensors longitudinally deployed in Elkhorn Slough, 
Monterey Bay, CA, to describe tidal wave evolution and distortion. Elkhorn Slough is 
short, reflective ('-'100%) estuary consisting of a narrow, gently sloping main channel and 
vast marsh and mud flats near the landward boundary. The amplitude reflection for the 
astronomical tidal constituent is -'90%, whereas the distortion, described by summing all 
non-astronomical tidal amplitudes, reflection is large (>125%) for all stations stating that 
the seaward tidal wave is more distorted than the shoreward tidal wave. It was found that 
the reflective time, defined as the time it takes the tidal wave to propagate to the 
landward boundary and back, varies as function of tidal elevation. Surprisingly, the 
reflective time increases with increasing tidal elevation contrary to the hypothesis that the 
shallow-water tidal wave phase speed increases with increasing tidal elevation, resulting 
in a reduced time. The landward end of the slough has vast low-lying marshes and mud 
flats that are inundated during higher tidal elevations, causing the effective (average) 
wave depth to decrease with increasing tide elevation, which explains why the tidal wave 
phase speed decreases instead of increasing. The elevation-dependent reflective time 
caused by a time-varying wave phase speed modifies the shape of the seaward tidal wave 
relative to the shape of the shoreward tidal wave, which explains the evolution and 
distortion of a tidal wave in a reflective slough. 
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I. 


MOTIVATION 


As the focus of U.S. Navy operations has shifted from the blue water to the brown 
water to meet the changing threats, it is important that the war fighter has a vast 
understanding of the physical processes associated with all environments to obtain a 
competitive advantage against the adversary. To gain a greater understanding of the 
littoral environment, the Office of Naval Research (ONR) has funded several, large, 
cooperative, research projects to study inlets in the U.S.. Specifically, in May-June, 2012 
an experiment at the New River Inlet, known as the River and Inlet Dynamics experiment 
(RIVET) was conducted. RIVET 2 is currently scheduled for the Columbia River, in 
May, 2013. The New River Inlet is approximately 200 m wide with a very large back 
bay. While river is in the name of the inlet, there is no river forcing the inlet is actually 
tidally forced. The Columbia River is a different system in that it is approximately 2 km 
wide and has strong river forcing. A third system which is different from the systems 
described above is Elkhorn Slough which is a short (10 km), shallow (3 m), highly 
reflective system, where the tidal wave is a standing wave. Erom a short experiment 
conducted in the NFS field studies (OC4210) in 2009, differences were noticed in wave 
characteristics between the shoreward and seaward tidal wave. Due to a lack of statistical 
confidence in the measurements and the increased desire by the Navy to gain a better 
understanding of the physical environment of the littorals, a small field experiment was 
conducted to explore the tidal wave characteristics of a highly reflective estuary. 
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II. INTRODUCTION 


As a tidal wave propagates from the open oeean through the inlet and upstream 
into the estuary, it distorts due to ehanges in water depth, channel geometry and frictional 
effects (Friedrichs and Aubrey 1988; van Rijn 2010). The tidal distortion is generally 
represented by overtides (tidal harmonics of astronomical constituents) and compound 
tides (nonlinear sums and differences of all tidal constituents) that the astronomical tides 
develop in shallow water (Dronkers 1964; Uncles 1981; Speer and Aubrey 1985; 
Friedrichs and Aubrey 1988; Nidzieko 2010). The tidal wave distortion is described by 
the tidal asymmetry, which is the inequality of the ebb and flood tidal currents. The 
inequality determines the net sediment transport and flushing patterns (Dronkers, 1986; 
Friedrichs and Aubrey, 1988). For example, an ebb-dominated system is characterized by 
a shorter duration, fast ebb current relative to a longer duration, slower flood current, 
which is often associated with ebb-tidal delta morphologic feature (Boon and Byrne 
1981). 

Most of the original work on tidal asymmetry has focused on estuarine systems on 
the U.S. east coast, where the degree of the distortion is primarily described between the 
interaction of the M2, semidiurnal lunar tidal constituent, and the M4, the first overtide of 
M2 (Friedrichs and Aubrey 1988). Nidzieko (2010) expanding upon the work of Hoitink 
et al. (2003) and Woodworth (2005) demonstrated that the tidal asymmetry for mixed- 
mainly, semidiurnal, tidal regimes (i.e., short estuaries on the U.S. west coast) are 
produced primarily through the combination of the astronomical lunar (Kl), lunisolar 
(01) diurnal tides and M2. The phase relationship between Kl, 01, and M2 induces an 
ebb-dominant tidal asymmetry because the lower-low water (LLW) follows the higher- 
high water (HHW) (Friedrichs 1995) and the phase of sea surface elevation and velocity 
are in quadrature (Nidzieko 2010). Though there is an ebb dominance at the estuary 
mouth, Nidzieko (2010) showed that the tidal asymmetry in short estuaries is modified as 
a function of distance upstream due to the increased amplitudes of the overtides and 
compound tides, consistent with the work by Friedrichs and Aubrey (1988). Nidzieko 
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(2010) provides examples of estuaries, where the upstream tidal asymmetry (distortion) 
beeomes flood dominated, and others where the ebb dominanee is enhaneed. 

An interesting question is raised about the behavior of the tidal distortion owing 
to the standing-wave nature of a short estuary. A standing-wave estuary implies that that 
shoreward propagating tidal wave is nearly perfectly reflected at the landward boundary. 
If the shoreward propagating wave becomes distorted as it travels landward, does this 
distortion decrease (or increase) as it reflected back toward the ocean? To examine the 
shoreward and seaward standing-wave structure of the tidal wave requires observations 
of both the potential (pressure) and kinetic (velocity) energy. However, most experiments 
tend to use pressure measurements solely, because pressure sensors are robust and 
inexpensive. 

To address tidal evolution and distortion, a field experiment was performed in fall 
of 2011 in Elkhorn Slough (ES), CA, which is one of the short estuaries discussed by 
Nidzieko (2010). Owing to the standing-wave structure, an array of four co-located 
pressure and velocity sensors were deployed throughout the slough, such that shoreward 
and seaward tidal signal can be separated. The separated signals allow for an accurate 
description of the shoreward and seaward tidal propagation and corresponding distortion. 
To date, this is the first time that this approach is applied to tidal wave motions using 
field observations. The amount of tidal wave reflection for the astronomical constituents, 
overtides, and compound tides are described. The characteristics of the shoreward and 
seaward tidal wave distortion and the mechanisms responsible for the distortion are 
discussed. 
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III. EXPERIMENT SETUP AND DATA COLLECTION 


ES is located at the center of Monterey Bay (MB) and is adjaeent to the landward 
edge of the Monterey Submarine Canyon. The Monterey Bay Sea Floor Mapping Lab 
(MBSFML) at California State University Monterey Bay (CSUMB) performed high- 
resolution bathymetric surveys of ES between September-Deeember, 2011 (Figure 1). ES 
eonsists of approximately a 12 km main ehannel that is 200 m wide near the mouth and 
narrows upstream to 90 m at Parson Slough (Broenkow and Breaker 2005). The length of 
ES is only 3.5% of the semi-diurnal tidal wavelength. The ehannel from the edge of the 
man-made jetties to the HWl bridge is approximately 4.5 m deep. East of the HWl 
bridge, there is an abrupt reduction in the water depth to 2.5 m (Figure 1). There is no 
ehannel dredging east of the bridge. There are a number of depth undulations along the 
ehannel, but the average slope is 1:4000. The deepest part of the channel is 10 m and is 
assoeiated with erosion by seeondary eireulation around Seal Bend (Figure 1). ES is 
surrounded by extensive Salieornia marsh and mud flats with a number of small 
tributaries (Broenkow and Breaker 2005). At low tide, the slough is a narrow ehannel 
with a large area of exposed dry marsh and mud flats, while at high tide the marsh and 
mud flats are inundated. 
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Figure 1. a) Google Earth image of Elkhorn Slough showing location of the 
4 ADCP locations with white dots. The white dashed line highlights 
the low-lying marsh and mud flats at the shoreward end of the slough, 
b) 2011 bathymetry map of ES provided by CSUMB SEME. ADCP locations 
are shown with black dots. Depth color-scale provided on the right. 


The slough is well-mixed, laeking stratification with water properties closely 
resembling MB (Broenkow and Breaker 2005). The dominant current direction in ES is 
streamwise fiow, and secondary lateral circulations generated by curvature in the bends 
of the slough are expected to cause flows to depart from the streamwise orientation 
(Wong 1989; Broenkow and Breaker 2005; Breaker et al. 2008; Nidzieko et el. 2008). 
Research has found that the tidal elevation and velocity in ES of the diurnal and 
semidiurnal tidal wave are approximately in quadrature as the result of standing-wave 
characteristics (Wong 1989; Broenkow and Breaker 2005). The inefficient water 
exchange at high tide due to the extensive Salicornia marsh and mud fiats and numerous 
small tributaries throughout ES are believed to enhance longer flood currents and shorter 
more intense ebb currents (Wong 1989; Broenkow and Breaker 2005). The wetting and 
drying of the marsh and mud fiats are hypothesized as the mechanism for enhancing tidal 
asymmetry (Wong 1989; Broenkow and Breaker 2005; Nidzieko 2010). 

Four 2-MHz Acoustic Doppler Current Profilers (ADCPs) with co-located 
pressure sensors were deployed longitudinally throughout ES channel for 66 days in 
Sept-Oct 2011 (Figure 1). The upward-oriented ADCPs were mounted on Im-high Ocean 
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depth (m) 








Science Group sea spider tripods. The four sensor locations were strategically placed 
throughout the slough. Station 1 and 2 (STl, ST2) were located on opposing ends of Seal 
Bend 2.5 and 3.5 km from the mouth, respectively. Station 3 (ST3) was located upstream 
of the connection to Parson Slough 5.5 km from the mouth. Note that Parson Slough is an 
approximately 50 m wide opening that connects to a complex set of tidal tributaries and 
marsh (Broenkow and Breaker, 2005). Station 4 (ST4) was located farther upstream 
approximately 500 m from Kirby Park and 7 km from the mouth of ES. The reflection 
point in ES is located approximately 12 km from the mouth. 

Velocities were obtained at 50 cm bin intervals sampled at IHz and averaged over 
90s before being stored. Owing to the minimal vertical variation, horizontal velocities 
were depth-averaged and rotated to the streamwise orientation. The relative across- 
channel variability described as the ratio of the across-channel velocity standard 
deviation to the streamwise standard deviation is 8%, 13%, 7%, and 7 % for ST 1-4. The 
flow at STl, 3 and 4 is considered streamwise dominant, whereas ST2 has more 
variability its proximity to the secondary circulation associated with the large curvature 
by Seals Bend (Breaker et al. 2008; Nidzieko et al. 2008). 

The ocean tide sea surface elevation for the corresponding ADCP deployment 
times was obtained from the NOAA tidal station (#9413450) located 23.5km south of the 
entrance to Moss Eanding in 7.3 m of water. The station uses a pressure sensor sampling 
at 6 minutes. 
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IV. METHODOLOGY 


To describe the evolution and distortion of the shoreward (+) and the reflected 
seaward (-) propagating tidal wave, the sea surface elevation (pT) time series are 
constructed from collocated pressure head (p) and depth-averaged streamwise velocity 
(u) assuming shallow water, 


ri- = \ll 



) 


( 1 ) 


where h is average water depth and g is the gravitational acceleration (Guza et al. 1984). 
This approach successfully described the cross-shore energy fluxes of infragravity waves 
in the nearshore surf zone environment (Guza et al. 1984; Elgar and Guza 1985; List 
1992; Sheremet et al. 2002). Tidal constituent amplitudes (A) and phases (0) for 
astronomical, overtides, and compound tides are computed for the p, and r\± time series 
using the T TIDE harmonic analysis program (Pawlowicz et al. 2002). Eor a 66-day time 
series record, T TIDE analysis resolves 35 tidal constituents describing 98% of the 
signal, where 2% of the signal is considered non-tidal residual motion. Only p and u are 
used to separate the signal, the across-channel velocity (v) is not considered. Large 
variability in the across-channel velocity (v) can bias p± estimates. Owing to the large v 
at ST2 and potential to bias p± estimates, results for ST2 will be omitted unless the 
variable of interest is p. 

In an effort to parameterize the tidal evolution and distortion, the tidal amplitudes 
are summed over two different frequency (f) bands, where the frequency cut-off between 
the frequency bands is 0.1 cycles per hour (cph). T TIDE amplitudes are summed for 
f<0.1 cph, which represent the primary astronomical forcing (i.e., M2, Kl, 01, among 
others), referred to as low-frequency motions (subscript Ifm). Amplitudes summed for 
1^0.1 cph represent tidal harmonics, such overtides (e.g., M4, among others), that are 
primarily used to characterize the tidal distortion, and are referred to as high-frequency 
motions (subscript hfm). Note that the total (subscript tot) respresents the combined sum 

of both the low- and high-frequency motions. Summed tidal amplitudes for the low- 
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frequency, high-frequency, and total motions are estimated for p and ri± (discussed 
below) and will be dennoted with subscripts. 

Tidal amplitude reflection (R) is the ratio of the square root of the seaward energy 
flux divided by the shoreward energy flux. 


R = 



l(A-( 

y2pg(A^)^ c^b \ 



( 2 ) 


where p is the density of sea water, is the shallow water wave group speed, and h is 

the channel width (Sheremet et al. 2002). When R=l, the tidal wave is perfectly reflected 
in the estuary and there are no energy losses in the system. An R<1 implies that there is 
tidal wave energy losses and R>1 implies there is an additional source of energy. Low- 
frequency, high-frequency, total R are computed for ST 1, 3, and 4 (Table 1). 



STl 

STS 

ST4 

Rtot 

95 

102 

100 

Rlfm 

90 

94 

92 

Rhfm 

125 

224 

231 


Table 1. Reflectivity coefficients for the total, Ifm and hfm for STl, 

ST3 and ST4. 
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V. RESULTS 


The behavior of the tidal wave is deseribed by tidal constituent amplitude spectra 
for measured pressure signals at MB and ST4 in ES (Figure 2). Low-frequency spectral 
amplitudes associated with the astronomical tidal forcing represent most of the tidal 
variance and are similar at both stations. The high-frequency spectral amplitudes are zero 
at MB, but increase farther upstream in ES. The high-frequency spectral amplitudes are 
integer multiple harmonics (i.e., M4, M6) of the astronomical forcing (i.e., Kl, 01, M2) 
(Figure 2). The amplitude spectra are consistent with previous spectra measured ES by 
Broenkow and Breaker (2005) and Nidzieko (2010). 
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Figure 2. Tidal constituent amplitudes estimated T_TIDE at ST4 in ES from 
pressure (black lines) and the long-term NOAA tidal constituents at the MB 
station. Vertical black dashed line at f=0.1 cph represents the frequency cut-off 
that separates the low-frequency motions to the left and the high-frequency 

motions to the right. 


There is a 16% increase in Atot from MB to ST4 in ES (Figure 3a) estimated by 
linear regression. This trend can be explained by the analytical solution for a co¬ 
oscillating tide. Lamb (1932) provides solutions for various simple geometric inlet 
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configurations. Of the configurations, ES is best described as constant width (b), linear 
varying water depth (h) resulting in a solution defined as, 


7 = c 




COS((Jt + s) 


( 4 ) 


where Jo represents the zero-order, first-kind Bessel function. 



, I is the slough 


length, ho is the slough depth at the end of the channel, a is the radian tidal frequency, x 
is distance from the landward boundary, C is the tidal amplitude, 8 is the phase 
difference. Equation 3 applied to ES geometries results in an 11% increase in amplitude 
from the mouth to the head of the slough, consistent with the measurements (figure 3a). 
While there is a difference between the model and measurements, the model describes the 
general co-oscillating tidal wave characteristic well considering the over simplification of 
the slough (constant width). 
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Figure 3. a) Squares represent Amt between MB and ST4 as a function of 
along-channel distance. The dashed line is the best-fit line for Atot estimated by 
linear regression. Solid line represents the Lamb (1932) analytical solution for 
co-oscillating tide for ES. b) Ap,hfm (squares) between MB and ST4 as a function 
of along-channel distance. The dashed line is the best-fit line for Ap,hfm estimated 
by linear regression, c) Rmt (triangles), Rifm (circles), and Rhfm (squares) as a 
function of along-channel distance. 
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The Ap^hfm (estimated with pressure) steadily inereases from 0 to approximately 
14 em between MB and ST4 (Figure 3b). The inerease in the tidal wave non-linearity 
with along-channel distance associated with pressure signal is consistent with previous 
research (Speer and Aubrey 1985; among others). However, the measured pressure signal 
represents the potential energy portion of a tidal wave. A purely progressive wave can be 
described solely by pressure or velocity because they are in phase. For a standing-wave, 
both potential and kinetic energy must be considered to effectively characterize the wave 
(i.e.. Equation 1). 

The reflection coefficients, R (Equation 2), for each station for the total, low- 
frequency, and high-frequency related shoreward and seaward tidal signal provide 
insights to the bulk tidal wave characteristics and its evolution. The Rtot ranges from 
95%-102% (Eigure 3c, Table 1). The small differences in Rtot are related to the fact that 
the propagating shoreward and seaward tidal can behave different at the different 
stations. Owing to this, STl represents the overall tidal behavior in the slough, because it 
is located the farthest from the upstream boundary and describes the reflection over a 
longer distance, whereas ST3-4 represents the tidal behavior farther upstream. The Rtot 
results imply that farther upstream the tide is more reflected than downstream. Eor STl, 
5% of the total amplitude is dissipated presumably by bottom friction (Speer and Aubrey 
1985) or secondary circulation associated with Seal Bend (Broenkow and Breaker 2005; 
Breaker et al. 2008; Nidzieko et al. 2008). In addition, the marsh and mud flats of Parson 
Slough located between STl and ST3 may contribute to the reflection at STl and its 
difference measured at ST3-4. ST3 and ST4 are closer to the reflection point and the 
complexities of the slough are much less, as there are minimal tidal tributaries and mostly 
marsh and mud flats. The high Rtot at each station suggest ES is highly reflective, 
consistent with Broenkow and Breaker (2005), Breaker et al. (2008) and Nidzieko 
( 2010 ). 

Rifmat STl, 3 and 4 are 90%, 94% and 92%, while R hfm at STl, 3 and 4 are 125%, 
224% and 231%, respectively (Eigure 3c, Table 1). Rifm indicate approximately 10% 
energy loss compared to the gain in energy associated with Rhfm. Owing to the fact that 
Rtot represents contributions from Rifm and Rhfm, and that Rtot is -100%, then the A^i+^hfm of 
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the seaward wave must be more non-linear than the of the shoreward wave. The 

larger AT,+^hfm of the seaward wave results in a Rhfm greater than 100% and aeeounts for 
the 10% energy loss assoeiated with Rifm. The larger Rhfm value suggests energy is 
aetually transferred from the low-frequency tidal motions (i.e., astronomical constituents) 
to the high-frequency tidal motions (i.e., harmonics) (Parker 1991). These are new 
findings previously unavailable from only pressure measurements. 
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VI. DISCUSSION 


The previous seetion highlighted the differenees between the shoreward and 
seaward waves using bulk statisties, but the bulk statistics do not describe and explain 
why the waves are actually different. Many researchers believe that tidal distortion in 
estuaries with a main channel and vast low lying marsh and mud flats is due to inefficient 
water exchange around high tide as the marsh and mud flats become inundated 
(Nummedal and Humphries 1978; Boon and Byrne 1981; Speer and Aubrey 1985; 
Friedrichs and Aubrey 1988 Wong 1989; Broenkow and Breaker 2005; Breaker et al. 
2008; Nidzieko 2010). Unfortunately, field measurements in this experiment were not 
obtained at the landward boundary to describe the effects of the marsh, but the time 
signals of the shoreward and seaward tidal wave are available to provide insight on the 
behavior that occurs at the landward boundary, which is examined next. 

The characteristics of a measured 24-hr time signal for yearday 269 of the 
shoreward (Figure 4, solid lines) and seaward (Figure 4, dashed lines) tidal wave for STl, 
STS and ST4 are evaluated. The shoreward tidal wave first propagates by STl (Figure 4, 
blue line) and then to STS (Figure 4, green line) and onto ST4 (Figure 4, red line), which 
is shown by the temporal offset between the stations. The temporal lag between each 
station represents the propagation time that is dependent upon the distance between the 
stations (e.g., the largest distance is between STl and ST4 resulting in the largest 
temporal difference). The seaward time signals (Figure 4, solid lines) are similar to 
shoreward time signals (Figure 4, dashed lines), except in reverse. The seaward time 
signal represents the tidal wave that has propagated to the landward boundary and 
reflected downstream. Therefore, ST4 represents the first station that the tidal seaward 
wave encounters, and then the tidal wave propagates past STS and then STl. The time 
difference between the shoreward signal and the seaward signal represents the time the 
tidal takes to propagate to the landward boundary and back, which is referred to as the 
reflection time. The reflection time increases with increasing distance downstream from 
the landward boundary. The reflection time is largest for STl. 
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Figure 4. The shoreward (solid lines) and seaward (dashed lines) sea surface 
elevation measured for yearday 269 (shown in relative hours) for STl (blue), 

ST3 (green), and ST4 (red). 

The time differenee between the shoreward and seaward time signal steadily 
increases (Figure 4) from low tide to high tide implying that the reflection time is 
changing and appears to be a function of elevation (total depth). The reflection time 
decreases from high tide to low tide and the response differs slightly from station to 
station (Figure 4). The reflection time for high tide maxima and low tide minima are 
estimated for the 66-day record for STl, STS, and ST4, resulting in four estimates per 
day. A linear trend is found between the reflection time and the sea surface elevation, 
which increases with increasing elevation (Figure 5). The linear slopes are 0.83, 0.99, and 
0.87 for STl, STS, and ST4 suggesting about a 1 hr time lag per meter of elevation gain. 
Though a linear trend is observed, there is a variability indicated by r -values of 0.59, 
0.69, and 0.60 for STl, STS, and ST4. The large scatter at STl (Figure 5a) may be related 
to influence of Parson Slough or Seal Bend. STS and ST4 interestingly have reduced 
variability at low tide (Figure 5b,c) suggesting that the reflection time below MSL is 
controlled by the channel and above MSL controlled by the more complex nature of the 
mud flats and marshes. 
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Figure 5. The reflective time lag estimated from daily high tide and low tide 
maxima and minima for a) STl, h) ST3, and c) ST4. Dashed line represents 
the hest-flt line estimated with linear regression. The linear slopes are 0.83,0.99, 
and 0.87 and the r^-values are 0.59, 0.69, and 0.60 for STl, ST3, and ST4. 

These results are initially eounter-intuitive, beeause it is expeeted as the depth of 
the main ehannel inereases (deereases) with a rising (falling) tide, the phase speed, 

^Jg{h + rj) , of the tidal wave would inerease (deerease) [Friedriehs and Aubrey 1988; 
Wong 1989; Broenkow and Breaker 2005] and deerease the refleetion time whieh differs 
with the observations. The total water depth (h+rj) assumes that the ehannel has vertieal 
walls. While this is relatively true for some portions of the slough loeated downstream, at 
the upper end of the slough the ehannel is surrounded by shallow marsh and mud flats 
(Figure 1), whieh beeomes inundated at higher tides (Broenkow and Breaker 2005). Sinee 
the ehannel eannot be deseribed by simple reetangle, the average (effeetive) water depth 
is estimated by 



where is the eross-seetional area and b is the ehannel width. Therefore, increasing tidal 

elevation results in a decrease in average water depth as water spreads out over the 
marsh. The result is a decrease in phase speed and increased reflection time as observed 
(Figures 4 and 5). This is consistent with the results found for a progressive wave system 

by Friedrichs and Aubrey (1988), where the volume of intertidal storage divided by the 
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channel volume is large in relation to the tidal amplitude divided by water depth that 
caused the wave at high tide to propagate slower than at low tide. 

The difference in reflection time as function of elevation modifies the seaward 
tidal signal relative to the shoreward time signal, whieh explains why the seaward tidal 
wave becomes more distorted. If the reflection time is eonstant, then the shoreward and 
seaward waves would be the same. However, this cannot happen even in a reetangular 
channel because the phase speed would increase with increasing depth (Friedriehs and 
Aubrey 1988). This suggests that the seaward tidal wave will always be more distorted 
than the shoreward tidal wave in a reflective estuary. 
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VII. SUMMARY AND CONCLUSIONS 


Co-located measurements of pressure and depth-averaged veloeity measured 
longitudinally throughout a short slough, allowed the tidal signal to be separated into 
shoreward and seaward tidal signals, providing new insights into tidal wave 
eharaeteristies. The overall tidal wave reflection in Elkhorn Slough, CA, is high (-'100%). 
The low-frequeney seaward tidal signal, assoeiated with astronomieal tidal eonstituents, 
tends to lose 10% of its amplitude relatively to the shoreward tidal signal, whereas, the 
high-frequency seaward tidal signal, assoeiated with tidal harmonics, has a large (>125%) 
inerease in amplitude relative to the shoreward tidal signal. The inerease in amplitude of 
the high frequency seaward tidal signal suggests energy is transferred from the 
astronomieal tidal eonstituents to the harmonies and overtides eonsistent with the 
findings of Parker (1991). The high-frequency amplitudes are responsible for the tidal 
wave distortion resulting in the seaward tidal wave being more distorted than the 
shoreward tidal wave. The meehanism responsible for the seaward tidal wave distortion 
ean be explained by the shallow-water wave phase speed, whieh is depth dependent. As 
the tidal elevation inereases, the small ehannel near the landward boundary fills up 
quiekly and water flows onto the adjaeent shallow mud fiats and marshes redueing the 
effective depth of the channel. The deereases in the effective depth reduces the shallow- 
water wave phase speed inducing a reflective time lag, defined as the time the tidal wave 
takes to propagate upstream to the landward boundary and back. A near linear trend was 
found for the reflective time lag that inereases with inereasing elevation for all stations. 
This trend effeetively modifies the eharaeteristie of the seaward tidal signal relative to the 
shoreward tidal signal, whieh eauses the tidal wave distortion. In summary, the refleetive 
behavior of the slough is responsible for the inerease in tidal distortion for the seaward 
tidal wave. 
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